Entanglement generation in harmonic chains: tagging by squeezing 
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We address the problem of spring-like coupling between bosons in an open chain configuration 
where the counter-rotating terms are explicitly included. We show that fruitful insight can be gained 
by decomposing the time-evolution operator of this problem into a pattern of linear-optics elements. 
This allows us to provide a clear picture of the effects of the counter-rotating terms in the important 
problem of long-haul entanglement distribution. The analytic control over the variance matrix of 
the state of the bosonic register allows us to track the dynamics of the entanglement. This helps 
in designing a global addressing scheme, complemented by a proper initialization of the register, 
which quantitatively improves the entanglement between the extremal oscillators in the chain, thus 
providing a strategy for feasible long distance entanglement distribution. 

PACS numbers: 03.67.-a, 03.67.Hk, 03.67.Mn, 85.85.+j, 42.50.Vk 



I. INTRODUCTION 

The development of reliable strategies for quantum 
communication and information transfer has gained, in 
these recent years, an increasing importance in the quan- 
tum information processing (QIP) panorama. The reli- 
able implementation of a quantum channel for the ex- 
change and distribution of information is indeed central 
in many potential QIP applications Q, Q|. Intuitively, 
one could think about a scenario in which the quantum 
channel and the processing device are two different en- 
tities which have to be interfaced at the right time of 
a given protocol. This implies the ability of switching 
on and off the interfacing interaction with sufficient de- 
gree of accuracy and a reliable control at the single-qubit 
level, which are very demanding requirements in general. 
On the other hand, the idea of exploiting collective in- 
teractions of intrinsically multipartite systems, governed 
by external potentials which globally address the entire 
register, has very recently encountered the interest of the 
QIP community 3]. A global addressing scheme offers 
advantages in terms of controllability of the device and 
protection from the decoherence channels unavoidably 
opened by any sort of local external intervention. 

Inspired by the progresses performed in the design, 
coupling and management of bosonic nanostructures, 
which can behave quantum mechanically 0, important 
efforts have been produced in order to better understand 
the role that multipartite systems of coupled bosons have 
in the transfer and propagation of quantum informa- 
tion 0. The application of global addressing tech- 
niques to systems of continuous variable (CV) bosonic 
systems is appealing for a less demanding implementa- 
tion of CV quantum information processing. 

In this work, we re-consider the issue of entanglement 
generation in a chain of harmonic oscillators coupled 
through nearest-neighbor spring-like forces induced by 
an external potential which addressed the whole system. 
One of the points of interests in our analysis is the role 
played by counter-rotating terms (present in the inter- 



action Hamiltonian) in the entanglement generation pro- 
cess. This point, anticipated by the studies in Q, is ana- 
lyzed here by a change of perspective. Instead of solving 
by brute force the dynamical equations ruling the evo- 
lution of the bosonic register, we look for a formal de- 
composition of the time evolution operator in terms of 
linear optics elements, following the successful route ini- 
tiated in We believe that this alternative approach 
clarifies the entanglement dynamics within the register 
and provides a more transparent picture of the role of 
the counter-rotating terms in such a process. Entangle- 
ment is found to be always present if the counter-rotating 
terms are included in the interaction Hamiltonian. How- 
ever, we find the degree of bipartite entanglement be- 
tween the first and last oscillator to be very small (a fea- 
ture which is evident, despite it has not been stressed, in 
the analyses in In order to quantitatively improve 
the entanglement settled between the ends of an open 
chain, we design a strategy based on proper initialization 
of the register (performed by locally acting on the state of 
the extremal oscillators only) and global addressing, fol- 
lowing the same lines depicted in quantum state transfer 
protocols 0- We show how, physically, this improve- 
ment is possible because of the symmetry properties of 
the bosonic system. 

The reminder of the manuscript is organized as follows. 
In Section [n] we introduce the interaction model here at 
hand. We discuss the technical tools used in order to 
derive effective decompositions of the time-evolution op- 
erator into linear optics operations. Effective all-optical 
setups can be introduced, which provide a visual picture 
of the evolution of an ./V-element register and we give 
an explicit example for a simple case. In Section IIIII 
the entanglement generated in an open chain is quanti- 
fied by means of the corresponding equivalent decompo- 
sitions. We show that, as long as only the quantum corre- 
lations generated by the counter-rotating terms alone are 
considered, end-to-end entanglement in the chain is not 
favoured. Strong quantum correlations, which never dis- 
appear, are found between the first and the second oscil- 
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FIG. 1: (color online) Diagrammatic scheme of the coupling 
configuration in Eq. Q . We consider N harmonic oscillators, 
in an open linear chain, coupled through nearest-neighbor in- 
teractions. The squeezing operations performed onto oscil- 
lator 1 and N are part of the tagging by squeezing scheme 
suggested in Section llVl 



lator in the chain. On the other hand, the entanglement 
between the first and the last oscillator is always very 
weak. A transparent physical interpretation of the time 
delay with which entanglement appears in the first-last 
oscillators subsystem is possible through the analysis of 
the corresponding all-optical setup. Section llVl addresses 
a way to improve the results discussed in Section II I II 
By considering the physical system as a fictitious two- 
terminal quantum black box, we show that simple local 
pre-squeezing of the first and last element of the chan- 
nel allows us to obtain several interesting effects. The 
end-to-end degree of entanglement can be quantitatively 
improved and any other bipartite 1 — > j quantum correla- 
tion (j = 2, .., N—l) can be correspondingly suppressed. 



II. THE MODEL AND THE EFFECTIVE 
DECOMPOSITION 

We consider TV oscillators labelled by j € [1, N] and 
arranged in an open linear chain. The coupling between 
the oscillators is provided by a nearest-neighbor spring- 
like force settled by an external potential. By including 
the free dynamics of each harmonic oscillator, the corre- 
sponding Hamiltonian reads 
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with qj = (bj + b'-)/y/2 and pj — i(bj ~ bj)/y/2 the po- 
sition and momentum quadrature operators of the j-th 
oscillator respectively and bj (K) the corresponding an- 
nihilation (creation) operator [8j- The coupling rates k's 
are real and time- independent. A sketch of the interac- 
tion configuration is provided in Fig. ^ The form of the 
coupling terms deserves some comments as it is straight- 
forward to see that each nqjqj+i in Eq. l|T)l. expressed by 
means of the annihilation and creation operators, reads 



Kqjq j+ -i 



i = ^(bX +1 + bM +1+ ^L +1 +m +1 ). ( 2 ) 



Eq. (0) includes co-rotating terms (bjb^ +1 + h.c.) as well 
as counter-rotating (CR) terms (bjbj+i + h.c.) [8|. 



In this paper, we treat the CR terms on the same 
footage as the co-rotating ones, analyzing their rele- 
vance in entanglement generation processes in interacting 
bosonic systems. In order to analyze the time-evolution 
of the chain, we look for an effective decomposition of 
the propagator U(t) = e~ chaint in terms of linear- 
optics elements. We order the quadrature-operators as 
x = (<?i)<?2i -,Qn,Pi,P2, ■■,Pn) T and divide H chain as 
H cham = H p cham + H q chain , where H q cham {H p cham ) in- 
volves only the g-part (p-part) of x. In matrix form 
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This is a tridiagonal matrix whose dimension depends on 
the number of elements in the register. The formal diag- 
onalization of H q h ■ guides us in expressing the g-part of 
Eq. J[J in a picture defined by eigen-operators which are 
linear superpositions of q^s. On the other hand, H p hain 
is already diagonal in the x basis and its form is not 
changed by orthogonal transformations. Therefore, we 
discard this part of Eq. Q from our explicit analysis 
and will include it only when necessary. 

The simple form of Eq. © allows for an efficient di- 
agonalization 0, which helps us in identifying a proper 
pattern of coupling operations for the decomposition of 
U. In order to keep our analysis general, we will refer to 
the well-known beam-splitter (BS) operator BSjk(6) = 
cxp[i6(qjpk — PjQk)] of its reflectivity sin 2 9 as a coupler 
operator because this term can be used for both optical 
fields and mechanical oscillators. In the eigen-operator 
basis we write 
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where Ej's are the eigen- frequencies of Eq. J3J and 

Of = Y^,k=i a jk9k (j = 1,--,A0 are the correspond- 
ing eigen-operators, expressed as normalized super- 
positions of the gfc quadratures with coefficient af k . 

The set {Ef,Of} is parameterized by the dimen- 
sion N of the chain. As an explicit example, we 
consider the first non-trivial case represented by an 
open chain of N = 3 where we have {Ef, E\, E%} = 
{uj/2, (u + V2k)/2, (uj - \/2k)/2}. We introduce the ma- 
trix of coefficients a^ k 
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Thus, the spectrum of H chain is a ladder, symmetric with 
respect to the bare eigen- frequency to/2. This is a general 
result for an odd number of oscillators: by increasing the 
dimension of the register, the spectrum of H^j? 1 (I G Z) 
is never degenerate ans consists of I different pairs of 
frequencies symmetrically shifted with respect to u>/2. 

Coming back to our example, the structure of the 
eigen-operators (5| (j € [1,3]) turns out to be very infor- 
mative in the research for a set of operations which can 
be used, starting from Eq. J3J), in order to get the diag- 
onal form QJ. Indeed, Eq. JSJ), suggests that a 50 : 50 
coupler operator BS13 simplifies the structure of the cou- 
pling terms leaving us with oscillator 1 being decoupled 
from the dynamics of the rest of the register. Then, a 
50 : 50 B23 operation will complete the diagonalization 
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with (jj which are the new quadratures after BS's to be 
put in correspondence with OjS. The matching with 
Eq. I0J is evident. Thus, the explicit diagonalization pro- 
cedure of the q-part of the chain's Hamiltonian gives us 
information about the pattern of BS operations which 
have to be applied to the bare expression Eq. @ in or- 
der to get Eq. In Fig. [21 we provide the sequence 
of BS's to apply in order to diagonalize the interaction 
Hamiltonian for N — 3, 4 and 5 (panels (a), (b) and 
(c) respectively). Straight lines represent 50 : 50 BS's, 
while curved ones stand for unbalanced BS's. Some com- 
ments are in order. First, there is a striking difference 
between the even and odd number of oscillators. For the 
odd number case, up to N = 7, a single unbalanced BS 



balanced BS's between the pairs of oscillators (1,2) and 
(3,4). A second important difference between chains of 
opposite parity will be highlighted later on. We stress 
that there might be other, inequivalent ordering of cou- 
pler operations which diagonalize the interaction Hamil- 
tonian. The choices presented in Figs. [3 (a) and (c) allow 
us to write the interaction part of the chain Hamiltonian 
in the special diagonal form corresponding to the associ- 
ation between oscillators and eigen-frequencies shown in 
the bottom part of Fig. [3 There, the set {Ef} is written 

as {uj/2 7 £i , ..,£jf_i} with the subset {£f} arranged in 
increasing order of frequencies (with k > 0) [Icj . 

After the action of the collective coupler operator 
B C oU, which collects the pattern (for a given N) dis- 
cussed above, the total Hamiltonian of the chain reads 

H q chain = (^/2)[(6f ) 2 +£jLi(jf ) 2 ]+Ef=i 1 ^i v (of ?■ 

Here, {Pj^} is the new set of momentum quadrature 
operators determined by the application of the coupler 
operations to {pj}- By specializing the discussion to 
the odd number of oscillators, the next step in order 
to find the decomposition of U (t) is the introduction of 
proper operations (acting on the elements of the regis- 
ter, oscillator 1 exluded) which balance the differences 
between oj/2 and £f. In the case of even N, these op- 
erations would involve the entire set of oscillators, with- 
out exclusions, as the eigen-spectrum of H^ ain (I € Z) 
does not include the bare frequency lu/2. Conceptu- 
ally, this balancing is an important step as it would 
allow us to look at the register as a set of new non- 
interacting harmonic oscillators. It is immediate to rec- 
ognize that this is possible through the use of single- 
oscillator squeezing Sj(sj) = exp[|sj(6f Pf + Pf Of )], 

which realizes — > e~ Sj Oj, Pj 
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is requiredm order to diagonalize H q cham . On the other thus formally write sl oU Bl oU U(t)B coll S coU = R coU (t) = 

SjURjiM*))' where Rj(4>j) = e l(t >J [{6 ? )2+( - p i N)2] is the 



hand, the BS pattern for the even number case appears to 
be more complicated already for N = 4, which is the first 
non-trivial configuration with even N, involving two un- 
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phase-space representation of a rotation operator of its 
angle cj>j and S co u = ^f-^Sjisj). By inverting the above 
relation, we find the formal expression 

U(t) = B co uS co uR co u(t)S\ oll B\ oU . (7) 

Analogously to the case of squeezing, the formal collec- 
tive rotation involves all the oscillators but the one asso- 
ciated to the bare eigen-frequency ui/2, which again is a 
specific feature of the odd N case. Moreover, the rota- 
tions R,j((f>j(t)Ys contain the entire time dependence of 
U(t). The balancing induced by the squeezing operations 
imposes, in general, time-independent conditions as it re- 
lates the squeezing parameters Sj 's to the elements of the 
set {£?}■ On the contrary, the formal identification of 



FIG. 2: (color online) Pattern of coupler operators (or BS 
operators) for the diagonalization of H c hain for TV = 3, 4, 5 
oscillators (panels (a), (b) and (c), respectively). Straight 
(curved) lines denote 50 : 50 (unbalanced) BS's. The time- 
order of the BS operations is such that BSi precedes BSj if 
i < j. The bottom part of the figure shows the correspon- 
dences between oscillators and eigen-frequencies induced by 
the corresponding BS's pattern. 
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coil Scoii with R co ii imposes that the rota- 
■j's carry an explicit time dependence. As an 



tion angles 

example, we consider again N = 3, where we have that 
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The decomposition Eq. J7J is a central result of our study. 
It allows us to provide a clear physical picture of the 
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FIG. 3: (color online) Equivalent interferometric setup corre- 
sponding to the decomposition of the time-evolution operator 
for an open chain of N = 3 oscillators. BS stands for beam- 
splitter (i.e. a coupler), R indicates rotation and S squeezing. 
M stands for a mirror. 



dynamics occurring within the linear chain, without ex- 
plicitly solving the dynamical equations of motion of the 
oscillators p|. Indeed, once the explicit form of B co u 
is found, one can straighforwardly infer the evolution of 
the oscillators configuration simply by considering proper 
squeezing and rotations. This is equivalent to design- 
ing formal interferometric setups which could be used 
for proof-of-principle experiments where, at least for a 
few elements, the effects of CR terms could be simulated 
and observed. Motivated by these arguments, in Fig. |3| 
we show the equivalent interferometer for N = 3. By in- 
spection, we see that this equivalent configuration results 
in concatenated Mach-Zehnder interferometers where the 
oscillators involved are subject to different squeezing and 
rotation operations. In going from N — 3 to N = 5, the 
overall concatenated structure of the setup is preserved, 
with just more oscillators being involved. This is not the 
case for N = 4, whose equivalent all-optical setup turns 
out to be more complicated than the one correspond- 
ing to N = 5 for instance, with squeezing and rotations 
involving the entire register, as already stressed. The 
intrinsic difference between the even and odd N cases 
should now be more evident. The most relevant discrep- 
ancy is caused by the absence of the bare frequency ui/2 
from the spectrum of H q chain ■ 

The second issue which has to be discussed here in 
relation to Eq. J7J is the role of the CR terms. It is 
easy to be convinced that the version of H c h a in where 
the CR terms are excluded would have a contribution 
having the form K^2j =1 PjPj+i- Now the p-part of the 
chain Hamiltonian is also non-diagonal, with the same 
tridiagonal structure in the quadrature operators basis 
discussed in Eq. ©. Therefore, the orthogonal transfor- 
mation which diagonalizes H^ h ■ (and the correspond- 
ing pattern of BS's) can be used in order to reduce the 
p-part as well, getting the same set of eigen- frequencies. 
The corresponding eigen-operators are superpositions of 
just the pj quadrature operators with the same numerical 
coefficients otj k 's appearing in Oj's. This implies that, 
for the odd N case, after the application of B co ii, we 
end up with Bl M H chain B cM = {uj/2)[{dff + (Pf) 2 ] + 



Ef= 2 ^[(Of ? + (P?n where Pf = £f =1 af kPk . Ev- 
idently, no squeezing is required in this case as the q 
and p-parts of the Hamiltonian are already balanced by 
the diagonalization procedure. Thus, the corresponding 
time-evolution operation could be immediately reinter- 
preted as the tensorial product of formal rotation oper- 
ators, one for each oscillator, showing that in this case 
the interferometric configurations sketched above are still 
valid: we only need to remove the squeezing operations. 

It should be clear, up to this stage, that the exchange 
of any information encoded in the elements of the bosonic 
register occurs entirely by means of the effective collec- 
tive operation B co u. The remainder of the decomposi- 
tion we have found, indeed, involves single-element oper- 
ations which do not mutually mix the oscillators. Thus, 
by considering co-rotating terms only, we can see that 
the structure of B co u, for a given N, is unchanged. This 
observation paves the way to the following consideration: 
the removal of CR terms from H c h a in prevents the spon- 
taneous creation of excitations in the system. In terms 
of the equivalent all-optical setups, this means that by 
preparing the register in a classical initial state, no inter- 
oscillator entan glem ent has to be expected, in this case. 
Indeed, in ref. 11] it is shown that non-classicality at 
the inputs of a BS is a fundamental pre-requisite for 
the entanglement of its outputs. In presence of CR 
terms, the single-oscillator squeezing provides the nec- 
essary non-classicality for inter-oscillator entanglement. 
We will come back to this point later, when the entan- 
glement generation is quantitatively addressed. 

It is worth comparing the decomposition in Eq. Q 
with what has been found for a star-shaped bosonic con- 
figuration |(|. In an open linear chain, B co u induces 
multi-body interactions between the element of the chain. 
In particular, from Fig. [3 we see that an exchange of in- 
formation is always required between the first and last 
element in a chain, a feature which holds regardless of 
N I n t ne star-shaped configuration, on the other 

hand, any exchange of information occurs via a preferen- 
tial way passing through the central component [|| . 



III. ENTANGLEMENT IN AN OPEN CHAIN: 
SYMMETRY OF THE VARIANCE MATRIX 

In order to investigate the dynamics of the entangle- 
ment generated among the oscillators in an open chain, 
we concentrate on Gaussian states and rely on the pow- 
erful tools provided by the variance matrix formalism. 
Indeed, the statistical properties of a Gaussian state, 
i.e. a state whose characteristic function is Gaussian, 
are entirely specified by the knowledge of its variance 
matrix. The variance matrix V is defined as V a p — 
({x a , xp\) {a, 0=1, ■■, 2iV), where x = q,p and, for con- 
venience, we have adopted the ordering of the quadrature 
operators x = (qi,pi, .., qN,PN) T - Throughout the pa- 
per, we assume that the Gaussian peak of each oscillator 
is at the origin of the respective phase-space. V is in one- 
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to-one correspondence with the characteristic function of 
a Gaussian CV state which, in turns, gives information 
about the state of the system 0] . 

When applied to an TV-oscillator input Gaussian state, 
the operations involved in Eq. JJJ give an output state 
which is also Gaussian. They can be formally described 
by means of the transformations 7r. (<^>j) = cos^jl + 
i sin <Ty for single-oscillator rotation and TgAsj) = 
e -sj<r z f or single-oscillator squeezing, where cr a (a = 
y, z) is the a-Pauli matrix. For two-oscillator BS we have 



such that 



H Jk (r 3 k,t jk ) = ( r J kl 



-r jk \ 
tjkl 



(9) 



where tjk, rjk stand for the transmittivity and reflec- 
tivity of the BS acting on elements j and k (with 
tfk + r jk = 1) E3- Explicitly, these transformations 
change an input variance matrix V to V Qj . = T^NT aj 
(a = S, R) for a single-oscillator 2x2 variance matrix 
and Vn , = V7r , for a two-oscillator 4x4 vari- 
ance matrix. From now on, we indicate with the final 
variance matrix resulting from the application of all the 
transformations involved in U(t) for a given TV. In this 
Section we focus the attention onto the case in which all 
the oscillators are prepared in vacuum state, so that the 
initial variance matrix of the joint state of the chain is 
V = ©j^lj. By using Eq. Q it can be shown that, for 
TV = 3, the final variance matrix reads 



'13 
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(10) 



where Li = 1 + C13 and L 2 = 1 + 2C13 account 
for the local properties of the oscillators while C12 = 
^Sj =2 (— l) J ' +1 Cj and Ci3 = \ Y?j=2 c ] describe the 
inter-oscillator correlations. We have introduced the el- 
ementary correlation matrices (which depend on the ef- 
fective squeezing and rotations of oscillators j = 2, 3) 



-e 2a * sin (<t>j) sinh(2Sj) I sin(20j) sinh(2sj 



I sin(20j) sinh(2sj 



2 1 

e 2sj gm 2 ( 



j) sinh(2sj) 
' (11) 

It is remarkable in Eq. (|10l) that the oscillators 1 and 3 
have the same local properties, which are different from 
those of the mediator oscillator 2. Moreover, the correla- 
tions between oscillators 1 and 2 appear to be the same 
as those between 2 and 3, which witnesses an evident de- 
gree of symmetry in the bosonic system ruled by Eq. • 
The proportionality of the correlation matrix C12 to the 
difference C3 — c 2 is important, in this analysis, and is 
in striking contrast with the inherent structure of the 
correlations between the end points of the chain. These 
observations will be crucial in the upcoming discussion 
relative to the improvement of the end-to-end entangle- 
ment. The structure of Eq. I|10|) is found to hold for larger 
registers. Indeed, as still manageable examples, we men- 
tion that for TV = 4 and 5 the decomposition of U(t) is 
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which extends the symmetry already manifested in Vj. 
In fact, the symmetry is a general property of Vy : it is 
straightforward to see that Vj exhibits symmetry with 
respect to the central element of the chain, whose local 
properties are unique in the system. The expressions of 
Cjfc's in terms of elementary correlation matrices analo- 
gous to Cj in Eq. (|11|) are, in general, quite cumbersome. 

We address the generation of quantum correlations 
among the elements of an TV-oscillator open chain as well 
as a simple strategy suitable for the improvement of the 
performances of this bosonic system as a long-haul entan- 
glement distributor. The Gaussian preserving nature of 
the transformations Tr.,- (</>j), Ts j (sj) and 7&- k (rji-,tjk) 
allows us to exploit the well-known necessary and suf- 
ficient conditions for the entanglement of two-body CV 
Gaussian states ^|>^J. The explicit object of our inves- 
tigation will be the evaluation of the bipartite entangle- 
ment between the first and the j-th oscillator in a chain 
of TV oscillators. Therefore, we consider the reduced vari- 
ance matrices v±j of the pairs (1, j) which are found from 
V/ by extracting the 4x4 submatrices (j — 2, .., TV) 
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V 2 j-l,2j-l 


V 2j~l,2j 


V v 2j>1 




V 2j,2j-1 





(13) 



As a measure of entanglement we use the logarithmic 
negativity which provides an upper bound to the en- 
tanglement of distillation [l^ and is strictly related to 
the extent to which a given state violates the Peres- 
Horodecki criterion for separability 16]. For bipartite 
Gaussian states, this entanglement measure can be eas- 
ily calculated starting from the symplectic spectrum of 
the partial transposition of the variance matrix v a b- In 
the phase-space, the partial transposition with respect 
to oscillator b corresponds to the time-reversal operation 
which flips the sign of the momenutm quadrature oper- 
ator of b. This can be represented by the action of the 
matrix P = 1 ffi cr z onto v a j,. We introduce the ma- 
trix Ytat, — S a © Sb, where E m = icr n , m (m = a, 6) is 
the symplectic matrix of oscillator m 14]. The symplec- 
tic eigen-values of v' ab — Pv a bP are the eigen-values of 
\iT, a bv' ab \, which are always equal in pairs. By calling j' n 
(ri = 1,2) the representative of each pair, the inequality 
min n (7^) > 1 is a necessary and sufficient condition for 
the separability of v a b- The logarithmic negativity A^ g 
is then evaluated as = Y^ n max (0> — ^°&2 7n) EH- 

The calculation of A^ for TV = 3, j € [2, 3] and 
k/u> = 0.1 leads to the plots shown in Fig.0J where the 
bipartite entanglement between the three oscillators in 
the open chain is plotted against the rescaled interaction 
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FIG. 4: (color online) Logarithmic negativity A„j (j = 2,3), 
for a three-element open chain, plotted against the dimen- 
sionless interaction time r = u)t, for k/ui = 0.1. The dotted 
line represents the behavior of A, 



A„g, the solid line A„ g . 



time r = cut. The choice for the ratio k/lo is dictated by 
the fact that, experimentally, a weak coupling regime of 
k < w is the only realistic situation 0, |g . The peri- 
odic behavior of the functions plotted is the signature of 
the time-dependence of the collective rotation R co ii(t). 
As seen in the symmetry of V^, we have A^ g = A^ 
(Fig. 01 dotted curve). The peak of A™ occurring at 
t = t* — 44.2 corresponds to |sin^| > 0.998. At this 
instant of time all the off-diagonal elements of Cj 's do not 
exceed ~5x 1CP 3 . On the other hand, it is evident that 
A^ 3 (Fig. 01 solid line) is smaller than A]^ g practically 
for any value of r, entirely disappearing at r*. Thus, 
despite the CR terms are responsible for the for free gen- 
eration of entanglement, a passive approach in which the 
bosonic register evolves freely without external interven- 
tion is evidently unsuitable for the creation of a reliable 
end-to-end entangled channel. On the contrary, almost 
all the quantum correlations within the system are local- 
ized among the nearest-neighbor oscillators (subsystems 
1 + 2 and 2 + 3). The trend is common to any other case 
we have checked: A]^ g can be almost an order of magni- 
tude larger than any other A^ (see Fig. for example) . 
Moreover, it is apparent that the behavior of each entan- 




glement function persists by enlarging the register. Only 
small modifications are observed in A^ g when N — > N+l, 
the most evident of which is that Ai^ becomes non-zero 
after an increasing time delay. This can be understood 
by considering the effective all-optical setup (as the one 
sketched in Fig.[3J): when all the oscillators are initially 
prepared in their vacuum state, the first set of BS's (on 
the left hand side of the figure) are ineffective as they 
superimpose |0) states. As soon as the squeezing of os- 
cillators 2, .., N is performed, the second set of BS's, on 
the right hand side of the figures, is responsible for the 
generation and propagation of quantum correlations. Ob- 
viously, the number of operations which precede the cou- 
pling between 1 and N, this latter carrying all the nec- 
essary non-classicality, increases with the dimension of 
the register, thus retarding the settlement of their en- 
tanglement. Again, the decomposition Eq. shades 
new light onto the important features of the entangle- 
ment dynamics throughout the system, complementing 
the results highlighted by previous analyses 



IV. TAGGING BY SQUEEZING 

For the purpose of creating an entangled state of the 
extremal oscillators in the channel, it is certainly desir- 
able to look for strategies which quantitatively improve 
the entanglement settled between 1 and N. Moreover, we 
would like to find out a way to make A^ dominant with 
respect to any other A^ g . The approach we are going to 
follow does not rely on local control over the elements of 
the chain between the first and the last oscillators. We 
assume that 1 and N are held by two spatially separated 
parties who can perfectly control the preparation of the 
respective oscillator and, if required, can measure their 
state. On the other hand, the interactions between the 
oscillators in the chain are set in a global way by a po- 
tential which collectively addresses all the elements at 
the same time. This approach is entirely within the rules 
of the global addressing strategies exemplified by quan- 
tum state transfer and phase-covariant cloning in quan- 
tum spin chains fA Il7| and by always- on computational 
schemes In this perspective, the chain is seen as a 
two-terminal device whose intermediate stage is embod- 
ied by the TV — 2 oscillators between the ending terminals 
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FIG. 5: (color online) (a): Logarithmic negativity A* J g (j = 
2,3,4), for a four-element open chain, plotted against the di- 
mensionless interaction time r = ujt, for k/oj — 0.1. The 
dot-dashed line is for A]^ g , the dotted line is for A^g and, fi- 
nally, the solid line represents A^ g . (b): Same as panel (a) 
but for N 



5. The thick dotted line is for A^g, the thin 

dot-dashed curve is for A*„ and the dashed one is for Alo- 
is 



Finally, the solid line shows A n 



FIG. 6: (color online) Scheme of principle for the end-to- 
end entanglement generation procedure. The two-terminal 
device embodies the bosonic quantum channel over which we 
assume to have no local control. On the other end, 1 and N 
are assumed to be within the grasp of two parties, which can 
arbitrarily prepare and measure the state of the respective 
oscillator. The interaction which determines the set of /t's 
globally addresses the A-element system. 
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1 and N. This central section is a black box whose dy- 
namics are out of the grasp. 

Intuitively, one would like to magnify the inherent dis- 
tinction of the pair of oscillators 1 and N, shown by 
Eqs. I|10|) and (|12|> . from the rest of the register. There- 
fore, any local action performed onto the ending termi- 
nals of the device in Fig. has to be designed so that 
the local properties of oscillators 1 and N still remain 
mutually equal. By considering the analysis performed 
in Section[n]and the role that non-classical states have in 
the entanglement by means of coupler operators \H\ , we 
look for an initial preparation of the register which can 
result in a quantitative increase of the end-to-end degree 
of entanglement. After a close inspection of the decom- 
position, we conjecture that single-oscillator squeezing 
operations onto 1 and N should improve the degree of 
entanglement between them. 

In order to demonstrate our conejcture, we address 
the case of N = 3 and we consider the preparation 
of an initial state whose variance matrix reads V = 
V sg ,i0l2ffiV sgi 3, where V S q,j = e~ 2rrT * d is the variance 
matrix of a squeezed state with its squeezing parameter r. 
The calculation of the logarithmic negativity for the sub- 
systems 1 + 2, 1 + 3 and 2 + 3 can proceed according to the 
recipe given in Section IlIII The corresponding degree of 



entanglement, from now on, will be indicated as A^ g ,taa 
(j,k 6 [1,7V]). The results, for r = 0.2, are shown in 
Fig. [7| (a) . Due to the symmetry of the particular initial 
preparation, it is easy to check that we get a final variance 
matrix having the same general structure as Eq. (|10fl . 
with suitably modified elementary correlation matrices. 
In particular, Fig. [7| (b) shows the time behavior of the 
matrix elements of Ci 2 = C23. Differently from what 
happens for an initially prepared vacuum state, when the 
ending elements are initially squeezed and for r ~ 46.8, 
all the elements {Cri)jk simultaneously become very close 
to zero (|(Ci 2 )i,i| ~ 6 x 1CT 3 , |(Ci 2 ) 2 , 2 | ~ 4 x 1CT 3 with 
|(Ci 2 )i, 2 | = |(Ci 2 )2,i| = 0), whereas C13 (at that value 
of t) becomes diagonal with matrix elements in the range 
of 0.1. This accounts for the improvement of the entan- 



glement settled between 1 and 3 with, correspondingly, 
A^ tag < 2 x 10~ 3 . The subsystem 1 + 3 has been tagged 
by the single-element pre-squeezing to be the preferen- 
tial pair of oscillators for the entanglement generation 
within the chain. It is worth stressing that this tagging 
procedure is possible in virtue of the symmetry exist- 
ing between the ending elements of the open chain. An 
analogous conclusion has been drawn in ref. |18| . where 
a totally symmetric A-body CV system has been con- 
sidered in order to point out the possibility of a unitary 
localization of the entanglement. In our case, however, 
the CV chain exhibits a degree of symmetry which is in- 
ferior to the one treated in ref. 18]. Different pairs of 
oscillators are characterized by different local and corre- 
lation properties, which makes the problem approached 
in this paper intrinsically different from the one in |18| . 
Nevertheless, we have shown that entanglement localiza- 
tion is possible with a lower degree of symmetry, which 
is per se an interesting point. 

We can generalize the choice for the initial variance 
matrix in the tagging procedure to the case of N oscil- 
lators by considering V = V sg4 © (©jl^lj) © "V S g,N- 
Again, an explicit calculation for the logarithmic nega- 
tivity can be performed, which leads to the plots shown 
in Figs. 0D (a) and (b), for the cases N = 4 and 5. In 
Fig. [S] (b) , the time-range has been restricted to the in- 
teresting region where A^ ta9 > A r \ J g ' taff , j € [2,4] in 
order to make the plot more transparent. It is evident 
that there is always at least one value of t at which the 
end-to-end entanglement dominates, making the tagging 
procedure effective. The amount of pre-required single- 
oscillator squeezing slightly depends on the dimension of 
the register and the plots in this paper show those values 
of r at which we have found a good trade-off between 
the degree of entanglement and the effectiveness of the 
tagging strategy. As we have stressed before, the even 
and odd cases are inherently different, as also witnessed 
by the fact that the value of r corresponding to an op- 
timized tagging operation is larger for N = 4 than for 





20 40 



FIG. 7: (color online) (a): Logarithmic negativity A^ tag 
(j = 2, 3) for a tagged three-oscillator open chain, plotted 
against the dimensionless interaction time r = cut, for k/u) = 
0.1 and r = 0.2. The dashed line is for A* 2 /" 19 = Af/ a9 
while the solid line shows A^' ta9 . (b): Time behavior of the 
elements of the correlation matrix C12 for the tagged chain 
of panel (a). 
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FIG. 8: (color online) (a): Logarithmic negativity A ] n ^ taB 
(j = 2, 3, 4) for a tagged four-oscillator chain against r = cut, 
for k/u = 0.1. The dot-dashed line is for A^/ a3 , the dotted 

14,iag 



line is for A^g' ta9 and, finally, the solid line represents A 
In this panel, r = 0.4. (b): Same as panel (a) but for N — 5. 
The thick black line is for A^'*" 9 , the thin dot-dashed curve 
is for A]^ tag and the dashed one is for K^ tag . Finally, the 
solid line shows A*g' t<19 . In this panel, r = 0.6. 
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N = 5. At the same time, A^' taff , for N = 4, is roughly 
proportional to r, while K]^' ta9 > r for all the odd N 
cases we have checked. 

V. REMARKS 

We have addressed the problem of long-haul entangle- 
ment creation in a register of bosons interacting via a 
global potential. The dynamics of entanglement can be 
clearly tracked via the effective decomposition of the time 
propagator in terms of simple linear optics elements as 
rotators, single-oscillator squeezers and couplers. This 
approach has allowed us to spot out a series of interest- 
ing features, characterizing the evolution of the quantum 
correlations settled among the elements of the register. 
As a result, we have been able to relate the conceptual 
role played by the CR terms in the entanglement gen- 
eration process to effective squeezing operations on the 



elements of the register. The usefulness of this analysis is 
also witnessed by the design of a tagging protocol for the 
improvement of the end-to-end entanglement and the si- 
multaneous reduction of any other 1 — ► j (j = 2, .., N— 1) 
quantum correlation in an chain of N elements. No lo- 
cal control over the central section is required: a proper 
preparation of the extremal oscillators and a collective 
interaction are sufficient to achieve the task. We believe 
this formal approach could be used in order to clarify 
other aspects related to the role played by the CR terms 
in entanglement creation, an issue which is certainly rel- 
evant especially in many problems of solid-state physics. 
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